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ABSTRACT 



The dynamics of phase locked loops have been studied by many 
investigators for the cases of either a pure sinusoidal signal or 
for a pseudo noise (PN) sequence. This dissertation presents an 
analysis of a phase locked system consisting of a sinusoidal carrier 
modulated with a PN sequence. The effect of the modulation index of 
the phase shift modulation and the design of a loop filter for coherent 
detection of a PN sequence are analysed. Additionally, the error 
probability for each digit, when a likelihood detector is used in 
front of the delay lock discriminator for a PN sequence, is also 
converted to the delay error of the phase lock loop. 
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TABLE OF SYMBOLS 



s(t)= an m-sequence with the number of digits M-2'’ -1 in a period, 
and the duration of a digit d, amplitude ±1 

0 = the modulation index of the phase shift modulation 

idt, = the angular frequency of the carrier 

A = amplitude of the incoming signal 

\j2Acos(i^c.t + 0s(t)) = incoming signal 

A 6 = expected value of the signal amplitude at the demodulator 
output 

G/»(**^)=the power spectral density of s(t) 
n(t)=* band limited white Gaussian noise 

n^(t)-an in-phase component of n(t) with respect to the carrier 
n^(t)»a quadrature component of n(t) with respect to the carrier 
N(<*^)*the Fourier transform of n(t) 

Y * one sided power spectral density of white Gaussian noise 
(watts/ cycle /second) 

F,(*o)-the loop filter of the demodulator 

Fa.(w)»the loop filter of the delay lock discriminator 

Be. = band-width parameter of F, (i^), (radian/sec) 

tdo = band-width parameter of Fa(«^), (radian/sec) 

(t)-- autocorrelation function of s(t) 

I^„('T‘)*crosscorrelation function of s(t) and n(t) 

ECxJ * the expected value of x 

6^*^- the variance of n(t) 
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CHAPTER I 



INTRODUCTION 



A correlation type detector or a matched filter is known to be 
optimum for detecting a signal in additive noise under a wide variety 
of conditions. The noise power is proportional to the effective band- 
width for the case of white Gaussian noise. We can Improve the S/N ratio 
by narrowing the effective bandwidth if it does not restrict the infor- 
mation contained in the signal. 

s(t)s(t+'t) 




Fig. 1-1 A Correlation Type Detector 



Consider the detector in Fig. 1-1. s(t) is the incoming signal and 
n(t) is additive white Gaussian noise. If we know the signal except for 
its phase, that is, the phase of the signal is the only Information 
contained in the signal, we can generate sCt+T) to multiply with the 
incoming signal. The low pass filter works like an averager in some sense 
as described later in section 4-3. and 
y(t/ , 

If s(t) and n(t) are assumed to be stationary. y(t,,T) becomes a function 
of onlyTior long T. and y(t) becomes approximately the sum of the auto- 
correlation function of s(t) "R^Cc)" and the crosscorrelation function 
of s(t) and n(t) "R^„(T)". 
y(T)» R^(tHr^„(t) 

If we can find a signal such that its autocorrelation function has 
a distinguished difference between the values of T=0 and XfO, we will 
be able to utilize it for phase detection. 
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A sine function ^ A cos cOtt has the autocorrelation function of 



I^(T) = A^cos and the phase locked loop for a sinusoidal function, 

like that shovm in Fig. 1-2, has been studied by many investigators. 

Among them, R. Jaffe and E. Rechtin showed an optimum filter design for 
a linearly approximated loop [ 5 j » A, J . Viterbi showed the transient 
response for 0 / having a second order polynomial in the nonlinear 
system, [lO] W. F« Roland defined the threshold criterion by the minimum 
of the cycles of slip and shovred its utility through experimental results 
in his thesis for the degree of M. S, [?] 

'{ZAcosiuJc.t t0,)+n(t) ACk„sin(6,-0^)t^Ck„n(t)sin(it^ct 



Low Pass 



' Csin(M4t +© 2 ) 




Filter " 




Voltage Controlled 
Oscillator 









Fig. 1-2 Phase Locked Loop for Sinusoid 
The binary signal has only two possibilities for each digit, ONE or 
ZERO, or -1 or -f-1, but sequences of binary digits show many interesting 
and ccmiplex characteristics. Some sequences represent codes for which 
maximum likelihood detection is possible. They are classified by the 
crosscorrelation function in each code set; orthogonal, biorthogonal and 
simplex. [2^ 

There are also some special group of binary sequences especially 
adapted for phase detection; the PN sequence, especially the maximum- 
length linear shift register sequence (m-sequence) , will be explained a 
little more precisely in Chapter 2, R, H. Barker defined "Group Synchro- 
nization" so as to pin-point effectively the origin in time for a binary 
sequence, and showed a matched filter, and also the characteristics and 
some patterns of this sequence, [l] 

The complementary binary sequences studied by S. Jauregui, Jr. used 
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the sane kind of matched filters and when used in a differential pair 
these eliminate small ambiguities in the output. (.6] Both sequences can 
be used in the form of discrete pulse sequences. 

PN sequences are adequate for the use of both continuously and 
periodically incoming signals. J. J. Spilker proposed a delay lock dis- 
criminator shown in Fig. 1-3. using the combination of two autocorrelation 
functions having two digit interval difference T. CsJ J. C. Springett 
showed a combination of clock and PN sequence to produce a controlling 
voltage proportional to the phase difference. f9] M. F. Easterling 
introduced a tracking system using a combination of clock and several codes 
and produced synchonization, code by code, separately. [ 3 ] Ih this system. 




Linear Feedback Shift 
Register 

Fig. 1-3 Delay Lock Discriminator for Binary Signal 



only the clock has a phase locked loop. The synchronization of the other 
codes is evaluated by a correlation indicator. For the system which needs 
long distance tracking and does not require quick reaction for the variation 
of the targets dynamics, this will be all right, but if we need automatic 
synchronization and very rapid acquisition, this will not be adequate. 

Since the second one has many weak ambiguous points between distin- 
guished phase lock points, we prefer Spilker 's delay lock discriminator 
rather than Springett 's. 

The periodic sequence is preferable for phase lock loops to the non- 
periodic because of the possibility of continuous control. 

This paper treats the phase lock system for a carrier phase shift 
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modulated by a PN sequence, as shovm in Fig, 1-4, s(t) is a maximum- 
length linear shift register sequence with the number of digits in a 
period M= 2*^-1, and the duration of a digit d, n(t) is white Gaussian 
noise and has the following narrow band approximation 
n(t)= n^(t)cos4t + n^(t)sini(^t 




Fig, 1-4 A Phase Lock System for a Signal Modulated 
by a PN Sequence 

Many papers have analysed each separate stage presented so far, but 
not the combined system. There are two problems to be solved for this 
system; what is the effect of 6 on the phase lock performance of this system, 
and if F,(w;) and F^(a)) have the same type of low pass filter as that derived 
by Jaffe and Rechtin, how can we determine adequate bandwidths for these 
filters? 0 determines the ratio of the I component and Q component of 
incoming signals, and also the necessary bandwidth for the input of this 
system. If 0 is large, the output modulating signal will be large, but due 
to the large phase error of the demodulating stage, the effective amplitude 
will decrease. The loop bandwidth of each stage is related to the tran- 
sient response for tracking and the required reduction in the noise output. 
This report will summarize those relations determined 1:^ other authors 
for each stage, and will further derive the threshold s/n ratio and the 
performance of this system. Chapter 2 considers the characteristics of 
m-sequences. Chapter 3 shov:s the demodulation of an incoming modulated 
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carrier and the output of this stage. 

The delay lock discriminator for an m-sequence will be analyzed in 
Chapter k with the effect of a limiter in front of this stage included 
in the performance. 

Chapter 5 summarises the system results and illustrates by an example 
the conclusion of the next chapter. This illustration will not show any 
particularly interesting value for 9 or other parameters. It is simply an 
example of the evaluation of system parameters and of design decision, 
but it is hoped that this will give a rough idea of the overall performance 
of such systems and will give some understanding of the role played by each 
parameter in the design of a system employing PN modulation of a sinusoidal 
carrier. 
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CHAPTER II 



MAXBin>I-LENGTH LINEAR SHIFT REGISTER SEQUENCE 

i'lio .•■'uLooo.rx’elation function of a binary sequence of length M, 
s(i)={a, a^|, vjhere a£ = ±l, is defined as 

M 

itj s modulo M 

If L(o)=i 

-&(j)=a4=l ( Kj<M) 

then s(i) is called a two-level sequence, and if, moreover, f^(j)=-pf 
for j^O (modulo M), the sequence is sometimes called a pseudo-noise 
(PN) sequence* Exactly, a PN sequence has to have the following three 
"randomness properties". [4J 

R-1 (the balance property) In each period of the sequence the number of 
one's differs from the number of ZERO'S by at most 1, 

R-2 (the i*un property) Among the runs of ONE'S and of ZERO'S in each 
period, one-half the runs of each kind are of length one, one fourth of 
each kind are of length two, one eighth are of length three, and so on 
as these fractions give meaningful numbers of runs. 

R-3 (the correlation property) If a period of the sequence is compared, 
term by tem, with any cyclic shift of itself, the number of agreements 
differs from the number of disagreements by at most 1. 

There are four types of PN sequences, corresponding to the length 

M. C2] 

1) M=2'^ -1 (ra-sequences) 

2) M=4t-1 is prime (quadratic residue or Legendre 

sequence) 

3) M =4t-l=4x*+27 is prime (Hail, sequence) 

4) M=p(p + 2) where both p, p + 2 are prime (twin prime sequence) 
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Among them the maximum-length linear shift register sequence (m-sequence) 
has the delay-and-add property and is generated by a linear shift register 
This is also convenient to generate in the phase lock loop, therefore in 
this report the n-sequence is used as a PN sequence. 




Fig. 2-1 n Stage Linear Shift Register 
In the n stage linear shift register shown in Fig. 2-1, the output 
sequence is always periodic with a period piM = 2"-l, and when period is 
M, that sequence is called a maximum-length linear shift register sequence 
This m-sequence has the delay-and-add property, that is, if the sequence 
is added term-by-term modulo 2 to a cyclic shift of itself, the resulting 
sequence is another cyclic shift of the original sequence. This property 
is possessed only by the binary sequence {a, .a^ , . . . ,a„} (ai=l or 0) which 
satisfies a linear recursion relation 



Pi 






for all k 
k-i : (modulo M) 

Ci* 1 or 0 

summation is modulo 2 

and M = 2" -1. Cz] m-sequence satisfies this. 

Consider an m-sequence s(i):s{a,, a^,...,a„} . The autocorrelation 

function of this sequence J^(k) is 

_^(k)=-p^^a{ea;,j^ 0, sum of modulo 2 

1*1 
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Using the del^-and-add property, 



= =-^{(# of ONE's)x(-l)+(# of ZERO' s)x(l)j 

= - ^ for kfO (modulo M) 

■^(0)=j5j2ai®a^ = 1 



i*/ 

(If we replace ONE and ZERO with -1 and +1, we can use multiple instead 
of © )o For the wave form s(t) which takes the value of +1 or -1 according 
to s(i) at time (i-l)d4t<id, the autocorrelation function J^(T) is, 

&(^ = ^ js(t)s(t+T)dt T«M 

ForT=0 or kd (k; integer) 

1 

^ Md I ^ 

S^s(t)s(t+kd)dt jdt g = -L 

If T=kd+(i' and 0<^i5^<d, then first (d-^) of each digit of s(t) 






corresponds with the last part (d-S') of s(t+kd), and the remaining S' of 
s(t) corresponds with the first part of s(t+(k + l)d). 

•^(kd +• ^ (Ai C^itk ^ I 

l_ 

~ M 

5a (T) is an even function and we can get the same result for any negative 
value of X d and the autocorrelation function of s(t) has the shape shown 
in Figo 2-2. 
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The power spectral density of s(t), G^(^) is derived by using the 
Wiener-Khintchine relation. 







When we consider the necessary bandwidth of the band pass filter for 
the modulated input \/2Acos(a;ct+0s(t)), where O<0<^ , we can analyae s(t) 
into frequency components. If we pick up one frequency, then that 
frequency component will be 




Since M is a large number, and the amplitude of each frequency 

component is supposed to be very small. A phase modulated carrier cOc. by 
a modulating frequency with small modulation index /?«!, is expressed 
as 

cos(K/tt+^sinWidi)« cos Wjt- ^sinu^^isinw^t 
This has frequency components of u>fj, The frequency components of 

the signal ^Acos(«<>tt+0s(t)) are derived in the same way as for small 
modulation index therefore in order to preserve almost all power 
of this incoming signal vre can pass only the frequency between 
because almost all the power of s(t) is included inside of 
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CHAPTER III 



PHASE LOCKED LOOP FOR CARRIER 
3.1 Coherent Demodulator 

Since s(t)=±l, i^e can separate the input signal -/^AcosCi^ct+^sCt)) 
into two parts. 

V^Acos (M^tt+^s (t ) )=/2Acos9 cosidjt+s (t ) \^As:lnOain UJJc 
If we can lock the phase of the local oscillator to th- c ricr frequency, 
v;e can get the modulating signal /2Asin0s(t). While V2Acos0cos t(^t can be 
used for that phase lock. Hence, there will be an optimum value of 0 to 
get the maximum signal to noise ratio in the output. 

Phase jitter in the local oscillator, caused by the output of 
\/2As(t)sin0, plus input noise causes the amplitude variation of signal 
output in the form of V2As(t)sin(0+ 9o) and acts to increase the noise 
input for the next stage. 

The crosscorrelation function of cos(uJtt+0s(t)) and cos«4t is 

f(X^=-:l=\ cos(«^ct+0s(t))cost4(tt-T)dt 

^ I J-j 

iC'^os6cos(2iUct -P Wcr)-s(t)sin0sin(2a^ct costOcXCOsS 

s (t ) sin0sin ftltX J dt 
= l-cos^cosiUtT-pj-sin^sina^cT 

— •|•cos0cos(^;T+ ^ ) 

<^= tan'*^tan0) 

<p is negligible vmen M is large and tan0 is not excessive and we can 
make a phase lock loop like Fig. 3-1. 

As far as this phase-lock loop is concerned, the first low pass 
filter can be neglected and letting the low pass filter be denoted as 
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B.v/. 




Fig. 3-1 Phase Lock Loop for Coherent Demodulation 
we can simplify this loop like Fig. 3-2. 



^Acos ( tL>c 1 1 X; + 0 s (t ) )+n (t ) ACkn,cos0sin (Xi~ \»)f ACk^s (t ) sinOcos (Xi - A o) 

® 



-^Csin(^yct+X.) 



-n(t)k/2Csin(wct -f-Xo) 



V.C.O. 



F,(«0) 



Fig. 3-2 Phase Lock Loop 

Noise input is assumed to be band limited white Gaussian noise, and 
using narrow band approximation 

n(t)=ni,(t)cosa>tt+' n^(t)sinidtt 

where n^(t) and n^(t) are normally distributed time functions varying 
slowly compared withtJc, and av. (nt )= av. (n^_)=0, the variances 
^=6Uj=N='^W W bandwidth (cycles/sec) 

>1 one sided power spectral density of 
white Gaussian noise (watts sec /cycle) 
VJhen we look at the vector of this input at the vicinity of the phase 
locked point in Fig. 3-3. \JZAsixid and n^ are found as the causes of 
phase deviation, and since they are independent, their effects will be 
evaluated in the sum of variances. 




Fig, - /cctcr Liagram for Phase Noise 
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E&r approximating 

sin( Ai~ Ao)= Xi - Ao 
cos(A^-X,)= 1 

in the vicinity of the phase locked point, we can linearize the loop as 
shown in Fig. 3-^» The output of the voltage controlled oscillator is 
the derivative of Ao and proportional to the input level. Thus ;^e can 
also modify them in Fig. 




Fig. 3“^ Linearized Phase Lock Loop 




Fig. 3-5 Linearized Phase Lock Loop 
Let k^Cgv=K 

Acos6=A| 

then we can get the loop of Fig. 3-5 • Tf we neglect the effect of 
0S(cO), the solution for optimum performance is shown in the paper of 
R. Jaffe and E. Rechtin. [5] 

Using that result, F,(u) is defined as F, (uf)= ^ 

the closed transfer function H(a)) is 



H(«>) = 



A,KF.M/ju> ^ B! + /z SoitJ 

f A,KFC^)/i^ + /a B. i w -I- ( j 



The network for F, («>) is approximately constructed by RC components. 
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Fron the linearized transfer function H(w), we can also derive the 



error 



6 fw)= 

= AiMf ' - ^ 






B* +/260/W 

For the input of frequency deviation 



(3-1) 



AU) 

B/ +\/2 6k/« - u)*- 



6w- 



= 0 
v2 



t> 0 
t <0 



Et - e-afdt - 

Max £(t) is derived by putting ^=0* 

tl = 

IrriAX^ Bo 



The effective bandwidth of H(tt») is 
= ^^pB(,= 1*06 Bo 



(3-2) 

(3-3) 

(3-^) 

(3-5) 

(3-6) 



(cycles/sec) 



(3-7) 



The phase error due to input noise is normal with zero mean and the 
variance is 

~ 21 

-00” 






dco 



3fl 



B = 



37R, 



2AWB 4^ 

3.2 Noise In the Demodulator Output 



(3-8) 



As shown in the previous section, the output near the phase locked 
point is 

Acos0sinQ>+As(t)sin0cos0o+n^(t) (3-9) 



var(n^(t)) = 

If we subtract the required signal As(t)sin0, then the remaining 
is unnecessary and is considered as noise, which is 

Acos0sin0o+ As (t) sin0(l~cos0o)-i-n^(t ) 
a) The distribution of 0© 

The phase error O# in eq.(3-9) is the sum of the contribution from 
the noise input and the self noise of s(t). 

Since H(w) is assumed to work in the linear region and is 
normal, the phase error produced by the noise input 0„ is approximately 
distributed normally. 

The distribution of phase error by the self noise 0s is derived 
as follows. 

Consider a randomness property of an m-sequence, "R-2" (the Run 
Property), "Among the runs of ONE'S and of ZERO'S in each period, the 
number of runs of length i is l/z'’ of the runs of each kind, and n = 
log^CM+l) is the maximtim length". VJhen we compare these distributions, 
we can see each digit has the probability of j to be +1 and to be -1, 
If we take a large number of digits m, then the probability it contains 
k of +1 is binomial distribution. 

P(k of +1 out of m digits) = 

Therefore if we take ra digits and count the difference of the number of 
-1 from the number of tl, the probability of the difference to be i is 

P(dif f erence - i)= | ( j ) 

(m and i should be both even, or odd, but since m is large, for a small 
number of i, we can adjust ra as m or m-i-1 according to the number of i 
for approximation). For a large number of ra, we can approximate it 
with normal distribution. 

Since a low pass filter is, in some sense, an averager, we can 



22 



conclude that Q, is also approximately normally distributed with zero 
mean and with variance equal to 0 Bod» (See Appendix l)/^ 

0^ and ^ have different kinds of input and they are independent, 
therefore, the probability distribution fxuiction of is also 

normal with zero mean and variance . 

b) The probability distribution of sin0o» and (l-cos 0,) where 4 has a 
normal distribution function. 

Considering the noise output by Qo, if we want to get good per- 
formance in the next stage , should be desired to be a very small value 
in the variance, if not, this system would not work satisfactorily. 

If 00 has a small variance, the contribution of 0oto Acos0sinG[) is 
approximately evaluated mainly by 9o neglecting the other terms of series 
expansion of sin0^, and for (1-cos 0„) lay Bi/z. 

Acos0sin0j= A^cos0 

As (t ) sin0(l-cos Oo)= ^As (t ) sin0 

0 * 

The first one has normal distribution and the second one is Chi- 
square distribution with decree of freedom 1. 

Looking at 0o which is the output of H(u)), we can see that it 
varies very slowly compared with each digit interval d, but since the 
ratio of maximum frequency deviation of the input to is the same as the 
ratio of the deviation of digit duration to d, and the filter parameter 
Bo is wider than the filter parameter in the next stage, 0<> varies more 
rapidly than the fluctuation of phase in the next stage. Therefore the 
effect of Oo to the next stage is considered as if Acos0sin6[> were added 
to n^ as random Gaussian noise in the approximation of (Acos0)0^, and since 
Asin0(l-cos 9o) is approximated by , at the output of the next stage they 
have the same effect as the decrease of amplitude by a factor like an 
average of many digits. 
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or 



The averaged power of As(t)sin0cos6j is 

A^sih*0 oos^de^^ e d9o 

•itt, V2n 

_ aW 0 rf .o , r I . 

^ [ / + ^ j;2^526^:c 

=A^sin^0[i^-ie-"‘^*/"J 

(l-cos0o)Asin0==-^ ■^* Asin0 



(3-13) 



h^) is distributed in %, distribution, and if v;e take sum of n 
V* 

digits, it becomes X„ 

E(^)= 1. var.(-^*)= 2/n 

And as n becomes large, the distribution of the average becomes very 
close to !! , 

av, (AsinO cos0o) = (1- -^°)Asin0 (3-1^) 

Comparing this vjith J , for snail values of 

+ ) =J 1- f' = 1- f'" 

The result shows good agreement wiili eq.(3-l4). 
c) Correlation of and n^t) 

From the argument so far, we could derive the approximate output of 
the disorl'inator , and it is sho>n as follov;s. 

'Hl+ e~^^) (Asin0)3(t> as the signal 

-f- (A cos0)0o + as noise (3-15) 

What is the distribution of noise as a whole? 

00 = + 9^ 

is uncorrelated v;ith 0,^ and so with n^(t). 

Therefore the correlation of 9q and i^^('t) is the correlation of Si, and 
n^(t). 



I 



In a linear system of transfer function H(k^), 
V2AO.S0 



h(t) = \/2B^e'^*^ cosj? t 



\T2 



for t>0 ( h(t)=0 for t<0) 
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Q,(t) - [“ n(t- T)h('r)dr 
Correlation of 9n and n^(t) is 

n(t+A) 0„(t)dt 

= lim ^ (^n(t+-A) ( n(t-'t)h(r)drdt 

7 > 0 & A V <T ^0 

= liBiT(i f**h(T) ^^n(ttA)n(t-r)dtdT 
= fh(r)lim \^n(t +\ )n(t- T)dtdf 

= j^h(T)R„(\-'OdT = jJh(r)^ (^X-t)dr 

for X= 0 

Ti V variances of 0<>and n^t) are 
t>«. Ufz k^ccs^& ^ ^ 

T'T ( 3 - 18 ) 

The normalized correlation coefficient f is 



y = 

Letting 



<S“>^ 

0(,Acos0 = X 
ng^(t) = y 
E(x)= E(y) = 0 
6i(^- A^cos^Q 6 ^^ 



(3-19) 



er*= g; 



n, 



Changing variables 



x+y = u 
y= V 



-<0 ;a2 

~ ®^C“ iCQ'+zffrxSy+s-/) ] jToo 

’"ft. 



where z = 



Jk 6i(5^ Cl 



.JL 






I a 2(/-/*) 

■2irV< “ f ^ 5V* +2f 6i«y r^* 

( > 



p ) -^ / r ^ 

L~ -z ( «/*■ t -2f <r« sy + J 
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(5‘;=6iV2f(S>(S>-(-S-/ 

U = Ilo(t) 

The output is now simplified as 
J id-f- e'^*^^) (Asin0^(t)+ no(t) 
np(t) is normal with zero mean and variance Sjf. 

The power spectral density of this noise is not flat in the bandwidth 
of l/d, because has the approximate bandwidth of B<,. Let the power 
spectral density of np(t) be Gn(<b) 

-03 

If we divide this by G„(wi , then we can get the effective band- 
width The low pass filter used in the next stage has narrower band- 

width than B^, then we can use ("^) as the approximation of G„(<i>) because 
the value of G„(< 4 >) inside of that bandwidth mainly contributes to the 
output. 

Let the Fourier transform of n^(t) be then the Fourier transform 

of (Acos 0 ) 0 o is ^^(ujyl^AcosB) 

(acos 0)%w)|^^== (AcosPlj^ i 0 

G,(o/)| = + = 2.917 

w-o " ai«<? ‘ 
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CHAPTER IV 

DELAY LOCK DISCRIMINATOR FOR A BINARY SEQUENCE 



4,1 Tracking and Acquisition Performance of a Delay-lock Discriminator 

This has been solved by J. J» Spilker, [Sj and here we would like to 
follow the general results. This section explains the tracking and 
acquisition performance, and the next section shows the effect of noise 
on the performance. 

If we take the difference of the output^ from two crosscorrelators 
which differ two digit intervals in Fig, 4-1, and separate it into an 
averaged value A^kDa(^) and another fluctuating part (self -noise) 
Aokn^(t,&), these are derived as follows. 




Shift Register 



Fig, 4-1 Delay Lock Discriminator 
The autocorrelation function of an m-sequence is already shown in 
Fig, 2-1, and from that figure we can draw D#|(€) as this. 
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D^(6)= E[s(t+T)s(tfdi-T)-s(t+T)s(t-d-i-T) 

= jR^(6-d)-R^(f+d)| (4-1) 

The self -noise in the output of the multiplier network is 
n^(t, e )= s(t^-d4)s(t-»-T)-s(t-d+T)s(tH')-D^(€) (4-2) 

For f=kd, k: integer, using the delay-and-add property, 

a) lc=0, (modulo M) 

n^(t.kd)- 5(t^-d1-T)s(t^■T) -s(t-d-t-T)s(t+T) 

= s ( tf T+d ) s ( t-t-T) -s ( t+T-d+d ) s ( tfT-d ) 

=s ( t+-T+jd ) -s (t+-T-d t-jd ) 

=s(ttT^od)-s(t+Tt(j-l)d) (4-3) 

b) k=±l, (modulo M) 



n^(t ,kd)= s(t+d^-Tfd)s(tfT)-s(t-d+-Ttd)s(t+T)-D(^(td) 



= (l-s(ti-T+id)- 



.s(t-hT+id)-l + 



Mt-I 

M 

Mfl 

M 



=+js(t-t-Ttid)- (4-4) 

c) k=-other integer 

n^(t,kd)=s(t-rT+rd)-s(t+Ttqd) (4-5) 

Further discussions of this self -noise, i.e. its power spectral 
density and the output from the loop vjill be explained in the next section. 
A) The output of the crosscorrelator network is 



kAo[D^(e)fn^(t,6)fno(t)/AoJ 



This is filtered by a low pass filter gfFa(«^) and drives the voltage 
controlled clock generator. This relation is now shown in the following 
Fig, 4-2. Obviously the control part of the crosscorrelator output D^(6) 

A 

is a function of T-T=^^, and T is the integral of the controlled frequency. 
T = ^ kAogc gf Fz (t<;) [D^(e )f n^ {(o, e)f n ^(w) /Ao j 
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Fig. 4-2 Simplified Delay Lock Discriminator 



In the locked-on state, €<d. 



D4(e) = 
e = T-T 



Mi- I g 
M 'T 



Using the same kind of filter as used in Chapter 3. 

-i-TZuJoiiO 









go — kAeg. gf 

X _ (m»* +'/ 2 K/o]o> r_e_ . <^»/Ao. 1 

A L it J 

The closed loop function H(w) is 
ufuj)- “^o*~ + ^ J 

aw- + w/ 

4 ^ (HtO/M J 

The value of T/d without noise is 



X=H(w)x 



Letting 



TM _ 






sin^ t 



^ — Wo ^ j:r t- 



max 



^__4_ 

4T 2‘1‘iiOo 



at t = 



•n: 

2{zU}o 



(4-6) 



(4-7) 

(4-8) 

(4-9) 



(4-10) 



(4-11) 

(4-12) 

(4-13) 



The critical value of uj, for e3Q)ected maximum value of a is to be 



decided frc«n nonlinear analysis. If we limit this as 
max. <?/d ^ 1 

then a/wo ^ 2.19 (4-14) 
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(^-15) 






' 2^ Wo^ 



dt 



B) Acquisition transient 



If we take onlyy from the previous equation, 

Jz Mo i IV 

>g- + 

I + J 



Jz |-N J /V \ ^ 



^0 I to_ T / Va<- \2. 
do ( «^0 J 



Dj(e) 



Let 



M 

Mtl 

A 



U>o 



x= y=7^ y-x’=T^ 7’^et 

implies 
C>d(^) 7 ^f is denoted D(x) 

D(x)=f 2-x 



X 

-2-x 

0 



f .D'(x)=| 1 

-1 



14 x4 2 
-14x41 
-2 4-x <-l 
othervjise 
-1 4 X 4 1 
-24x4-1 
otherwise 



then (y^ ('+ ^)DC^) 
y + y/g = ’x + x/g + D(x)f V2D'(x)x 



where x = 



x=-^ 



(4-16) 



go = 



g 



or 1 ^ X 4 2 



(4-1?) 



etc. 



x + (y +v/2D‘ (x))x-J-D(x)= y +y/g 

The slope of phase trajectory )T is 

Y-_x._ii_ D(X)t(j2 [/(x)+ /g 
® "T ~ ~ £ 

If we plot this performance on the phase plane x and x, D‘ (x) and 
D(x) are the function of only x and we can find the dividing line on 
x-=-2, x=-l, x = l and x=2. Let these regions be called A, B, C, D, E 
as shown in Fig. 4-3. 



(4-18) 

(4-19) 
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B 



Fig. 4-3 Regions in Phase Trajectcry 



- 'i 

t 



(4-20) 



Region A and £ 

D'(x)=D(x)=-0 

•• I • •• I • 

x+|x=y-efy ^ 

x= + V 
r+ '/j 

IT does not contain x, but is defined only by x and the input. 
Thus, if x = y at the initial point, then 



r= 

X = - 

y = x 




x»y-f-Xo x^s initial position (4-21) 

If x^y at T-0, letting the initial values of x, x, y, y be x,, x„, 0, 
y^ , and substitute them into the general solution of this equation 
x = c,-hCj e ^4- y 
and we obtain 

x = y + x,+ g(x, -y»)(l-e^^) (4-22) 

Region B 



D' (x)= -1 
D(x)= -2-x 

yr~ '2-xt('/s.-Jz)>[-W^ 
*_ (2 + X) + + i 

r+ 



(4-23) 



For a constant velocity input, i.e. y»0, the point having a constant 
T*» k draws a line . 

(kt^ - 2)x,-(x + 2 + y/g) = 0 
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If g is large enough to approximate l/g=0, then upper part of the line 
combining (x=-l, x= -0,518), (x®=-2, x=0) converges. This line is 




for constant T=k 

(k-i/2)x - (x + 2) = 0 



- 






\fz-k 



I 






sTT-k 



k<0 



the slope of the line= k = =- ... 

I 



therefore 



k^ -{2 k -1 = 0 

k - - n/3- 



' '/a 

Region C 

D(x)= X 
D'(x)= 1 

f (4.24) 

i t_a 

^ " r + '/j + /F 
Region D 

D(x)= 2-x 

D' (x)= -1 

I (4.25) 

V - -■2-i-x + » 

^ " r ± -sTE 

The line which divides the plane into converging and diverging parts is 
derived by the same way as that in region B. 



J 
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Fig. 4-5 Phase Trajectory in Region D 



The result obtained by CDC l604 digital computer ^ ^ ao and 
several initial conditions is shovm in Fig. 4-6, (Computer programming, 
see Appendix II). 

From this figure x=2.35» x=0 at the initial condition is the limit 
to converge. It Implies that -~<2.35* 

Time response of position error and velocity error are shown in 
Fig. 4-7, 8. From these figures the locking time is approximately 

6/o)o^ 27r/a;e 

The phase trajectory, x vs.tand x vs.Tfor g-*10 are shown in Fig, 

4-9, 10, 11. From these figures we can find the position error of 0.2, 
and the limit of y is approximately 2,0, Lock time is also 6lu}f>^Zir ItOo , 
4,2 Noise and its Effect on the tracking 

In this section, I would first like to discuss the effect of self- 
noise; second the effect of the input Gaussian noise, and finally, the 
effect of the limiter before this stage, 
a) Self-noise 

Self -noise is a function of time and phase error. Continuing from 
the last equation (4-2), power spectral density for the time error ^ 
being an integer multiple of d is derived as follows; 

n^(t,kd), where k=0 

The autocorrelation function of this self -noise is the sum of three 
triangular waves as shown. Ely the Wiener-Khinchin relation the power 
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spectral density is the Fourieij 
transform of this wave. 



G 






(«.,e»o) 



«t o 

— 2G^ (t«>) (l-C0Sii)d) 







(‘^.^=id)= G^M (4=27) 

Since n^(tpd) is the same as signal, 
ki= other integer 

^ transforming the autocorrelation function 
G„^(wp6»kd)= 2G^(«>)[l=cos(r=q)u>d] (4=28) 

We are interested in the effect of this noise in the linearized region 

-d^e^d. 



For the extreme case 6=±dp the power of self=noise in the output of 
H(u)) is 

=if U ;+ M ( -TF^7r^)„4 ' ^ ^ 

(4-29) 



g-i_ / 1 2 f 



^-00 
“d^Wod 

For the least extreme case , e = d 

/2jni_ li / • ^ .1 

' / A ' T*1 \*/ \ 

.75/ ( 

■ HAu)» 

b) The mean=square delay error caused by the input noise 

= 3-08d"(i^r%- ('»-3C) 

4,3 The effect of input limiting 

the previous discussion we approximated the input noise as 

2. 

Gaussian with zero mean and the variance The signal is a square 

wave of amplitude ±A<>« 
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Given a digit which has equal probability of +A or -A with additive 
Gaussian noise, and if we decide the signal is ±A according to the 
polarity of input, the error probability for each digit is derived as 
follows ! 

IVhen we defined A=Ji(l+ e’ Asin0, this was derived as an 
averaged value for some interval, hence if we want to get the probability 
of error for each digit, we can not use A© as the amplitude of the signal. 

The probability of error Pg. is 

Pg= P(s(t)= +l)P(output is negative/s(t)= +1) 

P(s(t)= -l)P(output is positive/s(t)= -1) (^-31) 

The conditional probability of each term is the same value, because 
additive noise has zero mean and the signal amplitude is the same in 
absolute value. 

Pg= 2P(s(t) = +l)P(output is negative/s(t)= + l) (4-32) 

The output of the first discriminator is 
Acos0sin0o+ As(t)sin0cos6.+ n|,(t) 

9o and Yi^ are Gaussian, with zero mean and 




they are not mutually independent, but the normalized correlation 
co-efficient fis 




From eq. (4-32) 

P(output is negative/s(t) =+l) 
=P [Asin(0t6,)-h n^< 0 j 



=P[n^^ -Asin(0^- )] 
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/ \ j r * 

where erf(x)=^)^-^ 

Consider a pulse train of ra digits, where each digit has probability 
of error Pg, then if all pulses are 1, the average of this pulse ti-ain is 
derived as follows. 

The number of error k in m digits distributes binormally, and if in 
is large enough and mPg is not small, then this distribution can be 
approximated by a normal distribution with mean= mPg , variance 
mPpCl-Pg.), In this case each error causes the decrease of 2/m in the 
average, and this decrease is normally distributed. 



The decrease of the averaged voltage v is, thus 



p(v) = 






Pe) 



(4-34) 



If the pulse train is composed with the same number of +-1 and -1, 
i.e. m/2 of positive pulses and ra/2 of negative pulses, then an error 
which occurred in positive pulses causes negative error 2/m in the 
average voltage and the same for an error in negative pulses, 

xs positive error voltage caused ty errors in negative pulses 
y: negative error voltage caused by errors in positive pulses 
Since the error for each digit is independent, x and y are independent 
because they are the error caused by only negative or positive pulses. 

p(y)= ! ) 

Since they are independent, the probability distribution of u = (x+y) is 



p(u)=[^f>^(X')p Cu-X^dX 






2n2F^0-f’^y^) iu-p,t-ru^p^-:cf)]4x 



2TT 6*2 







1 
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(4-35) 



This variance is equal to that in the positive pulse train. 

The effects these digit errors cause in the next binary discriminator, 
will be discussed in two phases, first at the locked point, and this can 
be considered as a phase jitter in the phase-locked position, and second 
at the limit where the phase error becomes ±d. This may be considered as 
the modification of curve. These two effects will show the analysis 

of the system. 

Vi/hen one digit error occurs in the input, how does this appear in 
the outputs of two crosscorrelators? If the system is locked in phase, 
^=*0, the outputs are also a cyclic permutation of that m-sequence and 
their phases differ by 1 digit. Since the autocorrelation function of 
this sequence is nearly zero for e»ld, the errors introduced in the 
outputs of the crosscorrelators do not have any relation to their polarity 
except the fact that they occurred at the same time. About the polarity, 
from the characteristics of the randomness of m-sequence, the probability 
that both errors are positive or negative is ^ each, and the probability 
the sum is zero is j* This implies the error probability is decreased 
by a factor of 2, and we can use P^/2 instead of Pg., and the change in 
the average voltage output caused by one error becomes double. 

Thus, we can modify the equation of the noise voltage as follows. 

^ ^ ~ ^^2if J ti(P*/2 . ) ( /- Ps/2.'> ~ f Z I (^“36) 

For the case of ^=d, one error always appears as negative, and 

another as positive of negative. This time we do not need to modify 
Pg, but we can add two effects. 

, , / , _ y _ I 

We have been discussing the effect of noise on the average through 
m digits so far, and now we have to make sure of the meaning of the 
average. 
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The inverse Fourier transform of H(^), h(t) is 
-^tr 

= cos^'t 

The response of this filter to a pulse of height 1 and duration d is, 
if we take the response after time d, 

^e‘%sin(^ t+^) 

Since the chance of error is very small, it will be improbable that the 
digit in error will be discovered. Therefore we can neglect the response 
of this filter in the digit duration and assume this as an impulse response. 

When we look backward for a time interval of T, standing at an arbi- 
trary time t|, if we use a real averaging device, the value at t, is, 
assuming there are k impulses to occur in T, 

«^(t-t^ )dt=^ . 

trj' 

If we use as an averager, the output average differs depending upon 

the positions of impulses, but when we think of the interval of each 
impulse, this is approximately distributed exponentially, and after taking 
many intervals T containing k impulses, we find k impulses are distributed 
in equal intervals as an average. 

Thus, the output average is the integration of that response. The time 
interval T should be determined for H(^), and looking at the response 
^e "/^sinC^t + , we can neglect the effect of the tail after 



m 3 ^ nJ a. 



(^-37) 



T — 

U^*^^ f e ^ Nf2sin( ^ t f )dt 



m = 



3lT\fz 



STTsTS: I 



srs rr 

V- M/« 






(4-38) 

(4-39) 



if 

Thus we can approximate the effect of H(»t>) by the effect of averaging. 
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CHAPTER V 



THRESHOLD FOR THIS SYSTEM 



5.1 Sunimai*y of the System Performance 

When we evaluate the performance of a phase lock system, maximum 
frequency deviation, or in other words, maximum target velocity that a 
tracking system can acquire, and transient time for phase lock, probabili- 
ty of cycles slipping and r.ra.s. phase error at the steady state will be 
considered. The problem of the cycles slipping is considered very im- 
portant for this system, because the acquisition time in A, E region in 
Fig, 4-3 can be very large. We can not allow cycle slip in the delay 
lock discriminator. Cycle slip is related to r.ra.s. phase error at the 
steady state. From the e3q>erimental results of V/, F, Roland, [71 we can 
find the linearity of B<> up to 20^-30 degree, and the expected 

o tr 

nximber of cycle slip per is 

6.10"^ for 6^. : 20® 

6.10'’'^ for 6-^, ; 25* 

2.35,10"^ for 6"^ ; 30“ 

If we use equation (3-3). the time for e =0.1 is ^ , 

This is for the initial position error of zero, and we can approximately 
use twice the value of it, , as an acquisition time frcsn the initi2il 
error of iz . 

Let the expected number of cycles slip per be p^ and the 
settling time be ^^/Bo . then as an average 

(5-1) 

is the number of cycles slip in one period. If one cycle of slip occurs, 

2TT 

in Md is in transient responce and it can cause longer transient 
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period of the next stage. 

For the cycles slip in the next stage, Spilker's paper shows that 
for delay error having 6i,^0.30d, a probability to lose the locked-on 
state in the absence of transient errors is negligible from experiments. 

While cycles slip determines the limitation for the values of the 
S/n ratio, d and the bandwidth parameter B<,, tc>o , the maximum frequency 
deviation restricts the lower limit of B<> and From the phase 
trajectory of the transient response of a phase lock loop for a sinusoid 
by A. ViterW. , the maximum frequency deviation at position error = 0 to 
be acquired is approximately 2.15Bo for this phase lock system concerned. 
If we restrict ^(t)l<l in the linearized system, 

^^2.19 from eq. (3-6) 

or -^<2.15 (5-2) 

Ho 

For the delay lock discriminator, from Fig, 4-6 and 4-10, 



This is 2. (5-3) 

Or, using the result of eq. (4-14), 

2.19 (5-4) 



The transient time for phase lock is decided mainly by B^ and , 
as follows, 

_2.tr 

‘*'lotk ~ phase lock loop 

+ s. 

Iftck ^t> 

These transient interval should be limited by the actual situation, 

5.2 An Example of the Calculation Procedure 

We are going to decide 9, B® and for a radar using this phase 
lock system. The specification of this radar is: 

Range $00 Miles 
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Carrier frequency 200 Megacycles 

Pulse duration, d 5 Microsecond ( 25 OO feet) 

Number of digits in one period M =^1200 

T— Md=5 Milisecond 



Maximum target velocity, y 3 I'lach, 3000 feet/sec 
Input S/N ratio less than 0 dB 

a) Lower limit of B^, 

^ = rooo 

Bt>^^ =4000 rad/sec 
for this B^, msec 

This will be short enough comparing with a period of Md—5 msec for 
the first acquisition and we do not need to increase B<». 

For the delay lock discriminator we put 



T(t)»adt 

^ 

1.4 2.0 



in eq, (4-11), 



0,6 rad/sec 

*l«ck=^ (5-5) 

for iJo=0,6 t, , s= 10 sec 

This is too large for this operation, and we need to make jJo large 
in the limit vdiere the cycle slip does not occur, 
b) The number of cycle slips per ^ 

( 5 - 6 ) 

B,= 4000 d = 5 . 10 ^^ M»1200 



Once a cycle slip occurs, it takes 1,6 msec to recover and this 
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becomes a loss of about one third of a period. Thus we want to limit 



-2. 



1^ in 3*10 . or r.m.s. phase error up to 25 • 

c) 

;„/R I I f >■“ ) 

^ 2A^c/^ 3^9 ' zce>a^e « ^ ^Tfxitoo k/o^/ 

where K is input S/N ratio 



S'>B.d0^ 



and ^^vs. d for various values of K are shown in Fig, 5-1 • 
d)4j = T(l^-e ®“^^)sin*^ 



A* 

Under the limitation of ^e^i0,19 ( = (25°)*^), we can plot ^ vs. 9 
for various K, as shown in Fig. 5-2. 

From Fig. 5-2, we can find 9 which maximizes with respect to 

K, and this is in Fig, 5-3. The maximized Af/k*' vs, K is shown in Fig. 5-^» 
e) The mean sequare delay error is the sum of the contribution of 
noise and self-noise. Taking the worst value of self -noise, from eq. (4- 
30) and (4-29), 

-^=3*08 jx(A^y^z^ 2- .2oo*/o-*7 (5-7) 

In Fig, 5-5, we showed vs, K with maximum value of At/ A*, for 

various values of /J., In order to avoid cycle slip, we limited ^/d^O.2, 
and the maximum value of «d<,vs K to satisfy this limitation, and t|^)^ vs, 

K for this oJoare finally shown in Fig. 5-6. 

PVom these figures, we can see that this system can not show good 
performance for the S/N ratio less than -15 dB. 

5*3 The effect of a limiter before the delay lock discriminator does not 
improve s/N ratio, because for low s/N ratio, ^approaches to ^ and phase 
lock control voltage D^(e-) effectively decreases to zero. 
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CHAPTER VI 



CONCLUSION 

The result in Chapter 5 shows that if we use the demodulator in the 
limit to avoid cycle slip, the amplitude of the output signal is mainly 
decided by sin 0 , and r.m.s. phase error 6^^ is also determined mainly by 
the noise input near the upper limit of 25 " 30 * • (5-7) shows that the 

variance of delay error O^^is the sura of input noise and self -noise and 
both are proportional to <^o, therefore for better s/N ratio the self -noise 
is dominant. We are interested more in the low s/N ratio region than in 
the high one , and we can also approximate 61. by the contribution of the 
input noise. 

When we want to design this phase lock system, we will find several 
opposing factors, 

l) The phase lock ability over the range of the frequency deviation 
limits the lower side of Bp, 

Z) The cycle slip bounds the upper limit of ^ and Bp, <Jo with respect to K, 

I'/ 

3) For a satisfactory acquisition time, we have to get lai’ge Bp and » 
The effect of d is related to the input bandwidth and d is also 

related to the range resolution of a radar. Small d will offer a good 
resolution, but on the other hand it costs S/N ratio and limits the input 
frequency deviation. 

We tried to solve the problem for a constant S/N ratio in Chapter 5. 
Actually input S/N ratio can not be constant in a system application, and 
we can estimate only the least desirable value or a range of s/N ratio. 
When we use this system for ranging, does not have much meaning in 
itself, but is very important as a cause of cycle slip. Looking through 
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these four factors, 2 and 3 are Especially significant in opposing each 
other, and we will finally decide Bo, and«^oby applying a weighting 
factor. The actual designing process must be decided according to each 
case, and we want to point out the fundamental relations between the 
above fotir factors, for they must be further investigated. 
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APPENDIX I 



SELF-NOISE 



The phase error output by signal s(t) itself is 0S('<>)H(‘*^) 
considering the relation shown in Fig. 3-5 • 
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Since M is very large, ^ is very small and we can approximate this 



relation by an integral form. 
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1 * 1 , -I^ -I 3 
Moreover, let g(2)» 

and consider contour integrals along the paths a, r, r, or n. 

Since B^d is positive, the contributions of the integral path H to the 
contour integrals jg(z)dz, and i Jg(z)e^ “ dz, and to j^g{z)~e dz 
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are all zero. The integral along 
the path T does not converge to zero 
for each integral as r-^-O, but the 
sum 



converges as follovfs. 



lira|yg(z)dz-T|g(z)e^^*^ dz-i j g(z)e"^^‘^ dz| 
»0 



Therefore 1 = 1, - 1^- Ij 



= 2 Trj\ 21 (residue for g(z)) - 21 (residue for g(z)e'* ) 

' UMP fiuD 



UHP 



UHp 



For I,, the sura of residues is -ffj, 

JJ_ 



+ i “S! (residue for g(z)e including z = 0) 

i~Hp 
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residue of z=0 is -jjjf- 
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( I- h£ 
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If Bod is small enough to neglect second power, 



I 



Thus 



= l + 1-3(1- -d- 

= 1.06Bed 
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APPENDIX II 

C014PUTER PROGRAM FOR EQUATION (4-18) 

PROGRAM ACQUIRE 

DIMENSION X(30),XI!OT(30),C(15) 

C(10)=l. 

1 CALL INTEGl(T,X,XDOT,C) 

IF (X(l)+2.0) 10.11,12 

12 IF (X(l)+1.0) 11,13,14 

14 IF (X(l)'l.O) 13,13,15 

15 IF (X(l)-2.0) 16,16,10 

10 DX-0.0 
DXDOT=0.0 
GO TO 2 

11 DX=-2.0-X(l) 

GO TO 3 

16 DX=2.0-X(1) 

3 DX.;OT=-:i.o 

GO TO 2 

13 DX=X(1) 

DXDOT=1.0 

2 ACCE=C(1) 

SPED=C(1)*T+C(2) 

ERROR=C (3 ) *SPEDfACCE- (C (3 ) +1 .414213»DXD0T) *X (2 )-DX 
XD0T(1)=X(2) 

XD0T(2)=ERR0R 
GO TO 1 
END 

60 



END 



For g=infinity 

time; from 0. to 9*0 » with interval 0,005 
C(1)=0. 

c(^;-y 

C(3)=l/g=0. 

ihitial value of X(l)“-2,0 
X(2)-y 

For g*=10 

time; from 0, to 1^,0, with interval 0.01 

C(1)=0 

C(2)=y 

C(3)=l/g=0.1 

Initial value of X(l)=-2.0 
X(2)=y 
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APPENDIX III 



SUBROUTINE INTEG 1 

This subroutine provides a simple but moderately flexible means 
for the solution of ordinary differential equations, employing integra- 
tion by the foxxrth-order Runge-Kutta numerical method. Input, and both 
print and graph output routines are built in. The step size of 
integration can be changed by appropriate data input and/or programming. 
Title: Runge-Kutta Solution of Ordinary Differential Equations 

with Built-in Input and Output Routines 
CO-OP ID: D2-NPS-INTEG1 (FORTRAN 60) 

Programmer; J. R. Ward 
Date : June 1964 
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